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ABSTRACT 



We report on the environmental dependence of galaxy properties at z = 2.15. We construct multi-band photometric data sets in 
the (proto-)cluster PKS 1138-26 field and in the GOODS field. We then fit spectral energy distributions of the galaxies with model 
templates generated with the latest stellar population synthesis code and derive physical properties of galaxies from the fits. To quantify 
the environmental dependence of galaxy properties, a special care is taken of systematic errors - we use data sets that have almost the 
same wavelength samplings, use the same code to fit SEDs with the same set of templates, and compare relative differences between 
the two samples. We find that the PKS 1138 galaxies have similar ages, shorter star formation time scales, lower star formation rates, 
and weaker dust extinction compared to the GOODS galaxies at z ~ 2. This trend is similar to that observed locally, suggesting that 
the environmental dependence of galaxy properties is already partly in place as early as z = 2.15. We show that the PKS 1138 galaxies 
assemble the bulk of their masses ~ 1 Gyr earlier than field galaxies, i.e., the galaxy formation depends on environment. Galaxy 
mergers should frequently occur during the first collapse of clusters and they might play an important role in driving the observed 
environmental dependence of galaxy properties at z = 2. 15. 

Key words. Galaxies : clusters : individual : PKS 1 138-26, Galaxies : formation. Galaxies : fundamental parameters 



1. Introduction 

The formation and evolution of galaxies in the Universe are 
dependent on environment in which galaxies live. In the local 
Universe, red early-type galaxies are the dominant population in 
rich galaxy clusters, while blue late-type galaxies are the domi- 
nant population in the low-density field. Not only galaxy proper- 
ties, but the formation epoch of galaxies also depends on en- 
vironment in the sense that cluster galaxies form earlier than 
field galaxies ( e.g..lKuntschner et al.l2002tlGebhardt et al.l2003l: 
iThomas et al.l l2005h . EarUer studies concentrated on nearby 
galaxies, but with the recent advent of large telescopes, environ- 
ment studies at z ~ 1 became possible. Interestingly, the envi- 
ronmental dependence of galaxy properties observed at z ~ 1 is 
already strong; clusters at z ~ 1 are dominated by red early- 
type galaxies (e.g ., Blakeslee et al. 2003; Nakataetal. 2005; 
iPostman et al."2005; Lidman et al. 2008; Mei et aI.II2009L but see 
also Cucciati et al. 2006) . Also, the formation epoch of clus- 
ter galaxies m easured at z ~ 1 is consistent with that ob- 
served locally (iGobat et al.l l2008h . Although ther e is a clear 
sign of galaxy evoltttion b etween z = 1 and (e.g. jElbaz et al.l 
120071; ICooper et al]|2008l) . one has to observe galaxies at even 
higher redshifts to fully quantify the environmental dependence 
of galaxy formation and evolution. 

Because of observational difficulties, only a few high red- 
shift clusters are known so far, the highest redshift cluster be- 
ing at z = 1-45 dStanford et al.l l2006). Higher redshift galaxies 
appear fainter and their rest-frame optical light migrates to the 
near-IR, where the sky background is brighter and it is challeng- 
ing to observe faint galaxies. Furthermore, high redshift clus- 
ters are poor clusters as they are still fast growing according to 
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the dark matter halo growth models (iPress & Schechtenll974l ; 
ISpringel et al.l 120051) . Such poor clusters are difficult to locate 
due to their weak density contrasts to the general field. There are 
a number of ways to find high redshift clusters, but one of the 
proven techniques is to look around high redshift radio galax- 
ies (Miley & De Breuck 2008). While not all the radio galax- 
ies are in over-density regions, man y of them host clear over- 
densities of galaxies around them ( Venemans et al.l f2007h and 
they are called proto-clusters. Although over-densities of red 
massive galaxies are not necessarily confirmed around them, 
they likely virialize and evolve to clusters at lower redshifts. 
Among the several proto-clusters reported so far, PKS 1 1 38-26 at 
z = 2.15 is one of the most promising proto-clusters for its clear 
over-density of sp ectroscopically con firmed galaxies by previ- 
ous studies dMilev & De Breuckl(2008l) . 

Early studies of the the PKS 1 1 38-26 ra dio galaxy at z = 2.15 
were performed bv iPentericci et al.] (Il997i il998 ). who reported 
a clumpy moiTjh ology of the radio g alaxy. Followed by these ini- 
tial observations. ' Kurk et al.l ( l2000h first reported an over-density 
of star forming galaxies around the ra dio galaxy. This region 
was then followed up by several authors. IPentericci et al.l (l2000h 
performed spectros copic follow-up o bservations of Lyman a 
emitters reported in iKurk et aP (l2000l) and confirmed 14 galax- 
ies close to the radio galaxy redshift. [Kurk et aP (|2004b'a) car- 
ried out further imaging and spectroscopic observations of the 
field targeting Ho- emitters and confirmed another 10 galaxies 
at the cl uster redshift. An X-r ay observation has also been per- 
formed ( Pentericci et al.ll2002h and follow-up spectr oscopy con- 
firmed at least 5 X-ray sources at the cluster redshift (ICroft et alJ 
2005). By now, there are more than 20 objects confirmed at 
the cluster redshift. Detailed analyses of the radio galaxy and 
the surrounding region with the superb resolution imaging with 
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HST have also been performe d (iMilev et alJ l2006t IZirm et aO 
l2008t iHatch et alJ l2008l l2009h . which added furth er hnes of 
eviden ce for the forming (proto-)cluster. In fact, IZirm et aTl 
('2008') r eported on t he forming cluster red sequence. Recently, 
[Pohertv et alj (l2009l) carried out a near-IR spectroscopic follow- 
up observation and confirmed two massive red galaxies at the 
cluster redshift. 

Given the the convincing over-density of galaxies and 
wealth of imaging and spectroscopic data available in the field, 
PKS1138 is an ideal sample to study the environmental depen- 
dence of galaxy properties at this high redshift. In this paper, we 
perform an extensive analysis of galaxies around PKS1138 to 
quantify the environmental dependence of galaxy evolution and 
formation at z = 2.15. 

The layout of the paper is as follows. We summarize our 
data in Section 2. We then describe details of our method of fit- 
ting spectral energy distributions of galaxies in Section 3. Before 
presenting our results, we perform extensive sanity checks in 
Section 4. Section 5 presents physical parameters of galaxies 
obtained from the fits as a function of environment at z ~ 2 
and Section 6 discusses implications of our results for galaxy 
formation. Finally, we summarize the paper in Section 7. 

Unless otherwise stated, we adopt Ho = 70km s"' Mpc"', 
Qm = 0.3, and Qa - 0.7. Magnitudes are on the AB system. We 
use the following abbreviations: IMF for initial mass function, 
SED for spectral energy distribution and SFR for star formation 
rate. 



2. Data 

We use two data sets for the purpose of this paper One is from 
the (proto-)cluster field PKS 11 38 to examine galaxies in an over- 
density region at z = 2.15. The other one is from GOODS for a 
field counterpart for comparison. We summarize the two sets of 
data in this section. 



2.1. PKS1138 

The field has been imaged by severa l instruments — [/-band 
with LRIS on Ke ck (Zirm etaD 12008), R band with FORS on 
VLT (iKurk et al. 



— ■■_ 112 000. 2004b), J, K, bands with MOIRCS on 

Subaru ('Kodama et al. 2007 ), 3.6-8. O^um with IRAC and 24;um 
with MIPS on Spitzer (Sev mour et al.l 2007). In addition to these 
published data, we used z-band data taken with FORS2 on VLT 
and //-band data with SOFI on NTT, which were obtained as a 
filler target of other programs and have not been published else- 
where. The z and //-band data were reduced in a standard man- 
ner. We further supplemented the data set with superb resolu - 
tion imaging by the ACS on-board the HST dMilevet al.ll2006h . 
We retrieved the pipeline reduced g and / band images from the 
Hubble Legacy Archive. Table[T]summarizes the data. 

The photometric zero points for the ground-based optical- 
nearlR images were o btained from the stan dard star observa- 
tions. As pointed out bv lDohertv et al.l ( 20091). the / an d K, band 
photometric zero-points used in iKodama et al.l (l2007h were off 
by ~0.3 mag. We re-measured the zero points from the stan- 
dard stars observed in the same nights and adopted those re- 
vised zero_2oints^or_the ACS images, we used the zero points 
from lSirianni et al] (l2005l). We have checke d all the zero points 
against stars from Gunn & Strvker (1983). We convolved the 
SEDs of stars from Gunn & Stryker ( 1983) with filter responses 
and atmosphere and derive synthesized magnitudes of stars. We 
then compared the observed sequence of stars with the synthe- 



sized sequence on color-color diagrams with various color com- 
binations. We found small zero point offsets with respect to the 
Gunn & Str yker (1983) stars and c orrected for them to better 
match with iGunn & Strvkerl (Il983h . The applied offsets were 
typically < 0.05 mag and the largest was 0.20 mag to the U- 
band. 

Objects were detected using sextractor ('Berli n & ArnoutsI 
1996) in the individual optical and near-IR bands because the 
seeing sizes vary from band to band and we did not smooth them 
to a common seeing. We used MAG_AUTO with aperture cor- 
rections assuming point sources for the optical and near-infrared 
photometry. For the IRAC images, we performed aperture pho- 
tometry using the /Tj-band image for object detections. We then 
applied aperture corrections to obtain the total fluxes. 

The aperture corrections were estimated as follows. We first 
constructed an average PSF image using bright, unsaturated stars 
in each band. Then we randomly distributed PSF objects in the 
images and repeated the object detection and photometry. The 
differences between the input magnitudes and measured magni- 
tudes were used for the aperture corrections. We also estimated 
detection limits for point sources and they are summarized in 
Table [T] The last column in the table shows the 50% detection 
limits. For images taken with ground-based facilities, apparent 
sizes of faint sources are small and magnitude limits for point 
sources are a reasonable proxy for extended sources. For the 
ACS data, the magnitude limits for point sources are extremely 
deep due to the superb resolution. We measured typical magni- 
tudes of the detected extended sources where their errors become 
0.3 magnitude (i.e., ~ 3cr limits) and adopted them as magnitude 
limits for extended sources. 

All the photometric errors were estimated from the sky noise 
in the same aperture sizes as used for objects. The photometric 
catalogs were cross-correlated with the /Tj-band catalog within 
the seeing FWHM of each band, and a /Tj-band selected catalog 
was produced. The final catalog contains UgRIzJHKs and IRAC 
3.5, 4.5, 5.8/iw photometry and the Galactic exti nction was cor- 
rected for in each band using the dust map from lSchlegel et al] 
(Il998i) . Stars are removed from the catalog based on their com- 
pactness and colors (we use a catalog in which stars are not re- 
moved only in Section 4. 1). 

Note that we have an IRAC E.Qi-im image, but we do not use 
it as strong dust emission (e.g., PAH emission at 1.6fj.m) from 
low redshift galaxies fall in this band. Our models described 
below do not include dust emission. We could in principle run 
photo-z without the 8/im photometry, select z ~ 2 galaxies, and 
then run the full SED fits to the selected galaxies including the 
8um photometry, which does not probe PAH at z ~ 2. But, that 
complicates the error analysis as it is not very straightforward to 
quantify how the errors propagate if we use two different sets of 
templates. 

We also have a MIPS image of the PKS 1138 field, but we 
only use it for a sanity check in Section 4. We apply a magnitude 
cut of Ks < 22.5 (~ lOcr) to the catalog to avoid any significant 
incompleteness effects. Also, we use only galaxies detected in 
more than 5 bands to ensure reliable SED fitting. 

2.2. GOODS 



We use data from the GOODS-MUSIC sample dGrazian et al.l 
2006; Sa ntini et all 2009) for a field counterpart of flie PKS 1 138 
data. There is a wealth of data in the GOODS field, but we re- 
strict ourselves to U35BVIZJHK, and IRAC 3.6, 4.5, 5.8/im data 
for this work so that we have almost the same wavelength sam- 
pling as the PKS 1138 data. This is a crucial point of this work 
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Table 1. Photometric data for PKSl 138. 



Band 


Instrument 


PSF size 


50% limit 


U 


LRIS 


1.2" 


26.9 


8 


ACS 


0.1" 


26.9 


R 


F0RS2 


1.0" 


25.3 


I 


ACS 


0.1" 


26.7 


z 


F0RS2 


0.9" 


24.1 


J 


MOIRCS 


0.7" 


24.9 


H 


SOFI 


0.8" 


22.3 


K. 


MOIRCS 


0.7" 


23.3 


'i.djim 


IRAC 


~ 2.0" 


23.1 


A.Sjim 


IRAC 


~ 2.0" 


22.5 


5.8yum 


IRAC 


~ 2.3" 


20.9 


8.0/jm 


IRAC 


~ 2.4" 


20.6 



- if we used more data in the GOODS field, then we would not 
be able to compare PKSl 138 and GOODS due to the different 
levels of systematic biases in the analyses. Note that we will use 
the MIPS photometry just for a sanity check in Section 4. The 
data covers a 78 arcmin^ field, which is limited by the //-band 
field coverage. We apply the magnitude cut of /T, < 22.5 to the 
GOODS catalog and use galaxies with detections in more than 5 
bands as done for PKSl 138. 

To make sure that GOODS is not a peculiar field, we de- 
rive galaxy density in the Subaru XMM-Newton Deep Field 
(SXDF) using the photometric data described in Finoguenov et 
al. (2009), which covered approximately 2,000 arcmin^. The 
density at z ~ 2 turns out to be similar to the GOODS den- 
sity; 0.57 + 0.02 arcmin"^ in SXDF and 0.64 + 0.09 arcmin"^ in 
GOODS. This suggests that GOODS samples a typical Universe 
at z ~ 2. 



- age, which is the time since the onset of star formation to the 
observed epoch. 

We allow r to vary between (single burst) and oo (constant 
SFR), and Ty between and 10. For the single burst model 
(t = 0), we assume no dust to mimic the passive evolution. For 
the other models, we do not assume any correlation between r 
and Ty. We adopt the two component dust extinction model of 
ICharlot & Falll (l2000l) . The actual amount of dust applied in the 
templates is not simply Ty but is dependent on star formation 
histories. We measure a V-band magnitude of a template and 
compare it with the V-band magnitude of dust-free template of 
the same age, t, and metallicity. The difference between the two 
is Ay. Note that we make fine model grids where z ~ 2 galax- 
ies populate. We implement e ffects of the intergalactic extinc- 
tion fol lowing iFurusawa et al.l (I2000.) . who used the recipe by 
Madaul ( Il995h and relaxed the allowed range of extinctions, as 
we explore the z > 2 Universe. Also, we use SWIRE AGN tem- 
plates from Polletta et al. ( 2007) as we mi ght expect an increasgd 
fraction of AGNs in t he PKS 1138 region (iPentericci et al.ll2002t 
iGalametz et"ani2009l) . 

We convolve the template spectra redshifted to various 
redshifts with response functions of the detectors and filters 
with/without atmosphere (we do not convolve atmosphere for 
data from the space) and generate synthesized fluxes. We use a 
logical constraint on redshift that the age of a template spectrum 
must be younger than the age of the Universe at that redshift. We 
have confirmed that we obtain essentially the same results if we 
do not use the constraint as we discuss below. The library of the 
synthesized fluxes for all the template is used to fit the observed 
SEDs of galaxies as detailed below. We use 2 million templates 
in total for both PKSl 138 and GOODS. 



3. The SED Fitting IVIetliod 

The idea of the work is to use the data from the two fields with 
almost the same wavelength coverage and sampling, feed the 
data to the same SED fitting code, fit them with the same set of 
templates, and discuss relative differences between the two sam- 
ples. We presented our data sets in the last section. We here de- 
scribe details of our SED fitting code. The SED fitting procedure 
follows the conventional x^ minimizing statistics. We prepare 
model templates and compare those templates with the observed 
data. We then analyze the fits and derive physical parameters out 
of the fits such as SFR and dust extinction. We start with prepar- 
ing templates and then describe details of the fitting procedure. 



3.1. Model Templates 

Th e model templates are gen erated using an updated version of 
the iBruzual & Chariot! (2003) population synthesis code, which 
takes into account the effects of thermally pulsating AGB stars. 
We adopt the Salpeter initial mass functioiQ (Salpeter 1955) to 
keep consistency with our previous study (iDohertv et al.ii2009i) 
and solar and sub-solar metallicities (Z - 0.02 and 0.008). The 
physical parameters that go into the templates are 

- star formation time scale (t) assuming the exponentially de- 
clining SFR (i.e.,SFR(t) cc exp(-f/T)) 

- optical depth of dust extinction in the V-band (ry) 



' The Chabrier initial mass function iChabrieJl2003h gives smaller 
SFRs and stellar masses by a factor of ~ 2 compared to the Salpeter 
IMF. Other than that, it gives the same results as those presented below. 



3.2. Fitting procedure 

We use the conventional ;t'^ minimizing statistics to fit the mod- 
els to the observed data. The fit itself is performed in fluxes, not 
in magnitudes. Non-detections are treated as zero flux with flux 
errors from the magnitude limits (here we take the 50% detec- 
tion limits). In this way, we avoid using upper limits in the fits 
and keep the x^ statistics simple. We impose a prior constraint 
that all galaxies have stellar masses lower than 6 x 1O"M0. This 
is a reasonable constr aint - the most massive galaxies at z ~ 2 
have ~ 5 X 1O"M0 (iKriek et al.ll2()08h . 

The best fitting template gives photometric redshift and 
physical parameters of galaxies. To be specific, we measure Zphot, 
T, Ay, age, SFR, and stellar mass from the fits. If the best fit- 
ting template is an AGN template, then we can not extract these 
parameters except for Zphoi because there is an non-negligible 
contribution of AGN to the overall SED. We only use their pho- 
tometric redshifts. 

An uncertainty on each parameter is estimated by marginal- 
izing over the parameter of interest and taking an interval of 
Ax^ = 1. However, all the parameters are correlated and this 
is just a rough error estimate. It is not easy to handle all the co- 
variances correctly, but fortunately, accurate random error esti- 
mates are not crucial for our purpose. It is because we discuss 
relative, systematic differences between the two data sets. We 
should have the same level of co-variances between the param- 
eters as we use the same code, templates and have almost the 
same wavelength sampling of SEDs. We therefore take this sim- 
ple error estimate, but we should bear in mind that our errors are 
not accurate. Fig.[T]shows an example of our SED fits. 
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Fig. 1. Sample SED fit. The points are the observed photome- 
try and the spectrum is the best-fitting template spectrum. For 
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4. Sanity Checks 

Before we present our results, we perform a few sanity checks to 
make sure that our results are robust. Here we discuss (1) pho- 
tometric zero points, (2) accuracy of our photometric redshifts, 
(3) accuracy of the SFR estimates from the SED fits, and (4) 
secondary burst models. 

4. 1 . Photometric Zero Points 

We cannot avoid random errors in the photometric zero points, 
but it is important to make sure that we do not have any system- 
atic zero point offsets as a function of wavelength, which can 
cause color stretches in the color space. For example, if the zero 
points are systematically fainter in the blue bands, we might get 
higher extinction and/or older ages from the SED fits. 

The easiest way to check the color stretch is to look at the 
sequence of stars in color-color diagrams. Fig. |2]is one such plot. 
The plot shows z - K, vs. u - z of the detected objects brighter 
than Ks = 22.5 in the two samples. One can easily identify the 
stellar sequence in the lower part of the plot. There is no strong 
systematic offset between the sequence of stars in the PKS1138 
field and that in the GOODS field. We have also checked various 
color combinations and found that the offsets are < 0.1 mags. 
To accommodate with any possible zero point eiTors, we add 0. 1 
mag in the quadrature to all the magnitude uncertainties in both 
PKS1138 and GOODS catalogs before performing the SED fits 
described in the last sectiorQ. 

4.2. Photometric Redshifts 

Another good check will be to look at the accuracy of photomet- 
ric redshifts. We want to make sure that we can cull the majority 
of the galaxies at z ~ 2 with a reasonably small amount of fore- 
/background contamination. 



^ The addition of random errors is not an ideal way to remedy the 
systematic errors. But, we apply the additional errors to make sure that 
the systematic errors do not dominate the overall error budget. The price 
we have to pay is that the additional errors slightly degrade the accuracy 
of photometric redshifts. 




Fig. 2. z- Ks plotted against u-z for objects brighter than ^5 = 
22.5. The open squares and dots show objects in PKS1138 and 
GOODS, respectively. Note the sequence of stars in the lower 
part of the plot. 




n PKS1138 - 
o GOODS 
X Y >3 

3 4 



spec 



Fig. 3. Zphot plotted against Zspec- The dotted line shows the 
Zspec = Zphot relation. The squares and circles are Ks < 22.5 
galaxies in PKS1138 and GOODS, respectively. The crosses 
show objects withal > 3 both for PKSl 138 and GOODS. 



Fig. |3] compares our photometric redshifts (Zphot) with 
spectroscopic redshifts (Zsnec)- The spectroscopic redshifts in 
PKS1138 a r e from Pentericci et al. (2000), Kurk et al. (2004^, 
ICroft et al.l (|2005), and Doherty et al. (2009) and those in 
GOODS are from the GOODS-MUSIC catalog dGrazian et all 
2006 and references therein). The correlation is reasonably good, 
although there are a number of galaxies with wrong photo-z. We 
find that galaxies with poor fits (the reduced ;t^^ larger than 3) of- 
ten have wrong photometric redshifts. We remove those galaxies 
from the main analysis. Our results remain unchanged if we do 
not apply this xl cut as we quantify later. 

A fraction of photo-z outliers (\Zspec-Zphot\/(i +Zspec) > 0.2) 
is 5% in GOODS over the entire redshift range. Our spectro- 
scopic sample for PKS 11 38 is strongly skewed at particular red- 
shifts, and we do not quote the numbers for comparison here. 
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The photo-z selection criterion for the main analysis is a trade 
off between the completeness and contamination. The selection 
adopted for the main analysis, Zphoi - 2.15 within 2cr, is pri- 
marily motivated by the rate of recovering the spectroscopically 
confirmed PKSl 138 members (it recovers 6 out of 7 confirmed 
PKS1138 members brighter than ^5 - 22.5). We have con- 
firmed that our results presented below do not change if we 
change the selection criterion within a reasonable range (e.g., 
\Zphot - 2.151 < 0.2). The contamination of fore -/background 
galaxies in our sample used in the main analysis is 36% in 
GOODS, which is estimated from the spectroscopically con- 
firmed galaxies at \zspec - 2.15| > 0.3 with photometric redshift 
consistent with Zphot - 2.15 within 2cr. Note that the numbers 
quoted here should be considered rough estimates because the 
spectroscopic redshifts are collected from various surveys and 
the sample is quite heterogeneous. As we will show later, con- 
tamination due to wrong photometric redshifts (e.g., foreground 
galaxies at Zspec ~ 1 but with Zphot ~ 2) do not strongly affect 
our conclusions. We note in passing that 6 out of 50 galaxies in 
PKSl 138 and 8 out of 50 galaxies in GOODS that are used in 
the main analysis are spectroscopic members. 

It is also important to verify that the parameters derived from 
the SED fits are not strongly affected by the uncertainties in the 
photometric redshifts. We take into account the random errors 
when we discuss differences between PKSl 138 and GOODS in 
the main analysis, but one may worry that photometric redshifts 
may introduce systematics in the SED fits and in the derived 
parameters. 

To quantify this, we fit SEDs at redshifts fixed at Zspec for 
the spectroscopic objects. Fig. |4] compares the parameters de- 
rived at Zspec and those with redshifts as a free parameter We 
observe a good correlation between them. In particular, there is 
no strong systematic offsets as shown in the top panel of Tabled 
where we show the median differences between the physical pa- 
rameters derived at Zphot and those at Zspec for the spectroscopic 
members at Zspec ~ 2.15. The values in the brackets show the 
rate at which the two parameters agree within Icr. Spectroscopic 
redshifts have narrow errors on all the parameters, but photomet- 
ric redshifts do not introduce any systematics. We will focus on 
systematic differences between the two samples in the next sec- 
tion and all the differences we observe are larger than the offsets 
listed in Table |2l Therefore, the uncertainties in the photometric 
redshifts do not affect our results. 

We will not use parameters derived at Zspec in the main analy- 
sis because we would like to show what parameters the contam- 
inant galaxies (i.e., objects with wrong photometric redshifts) 
typically have and to show that they do not affect our conclu- 
sions. Also, the spectroscopic samples are strongly biased to- 
wards particular types of galaxies and we do not want to put 
a weight on them. For example, many of the spectroscopically 
observed galaxies in P KSl 138 a re either emission line objects 
or X-ray sources (Pent ericci et a l. 2000; Kurk et al. 2004b) and 
they are not necessarily typical galaxies. We use photometric 
redshift as a free parameter for all the objects so as not to in- 
troduce any biases. We have confirmed that our results do not 
change if we use the parameters derived at Zspec where available, 
as expected from Fig.|4] 

Finally, we show that the logical age constraint - the age of 
a model template must be younger than the universe at a given 
redshift - does not change our results. The middle panel of Table 
|2] quantifies the effect of the age constraint on the derived pa- 
rameters. We observe no systematic differences and the derived 
parameters of most galaxies agree within Icr. We can therefore 
safely apply the constraint in the main analysis. The reason why 



Table 2. Differences in the physical parameters derived with 
other models or constraints than the fiducial ones. 



Parameter 




PKS1138 


GOODS 


age [Gyr] 




+0.13 (75%) 


0(81%) 


SF time scale 


[Gyr] 


0(100%) 


(70%) 


SFR [Mo yr- 


] 


-0.08 (88%) 


-0.6(81%) 


Av [mag] 




0(100%) 


(88%) 


age [Gyr] 




(92%) 


(96%) 


SF time scale 


[Gyr] 


(90%) 


(100%) 


SFR [Mo yr' 


] 


(76%) 


(92%) 


Av [mag] 




(80%) 


(97%) 


age [Gyr] 




(70%) 


(79%) 


SF time scale 


[Gyr] 


(72%) 


(80%) 


SFR [Mo yr- 


] 


-0.05 (56%) 


(73%) 


Av [mag] 




-0.50 (58%) 


-0.75 (37%) 



Notes. The top panel shows the median differences between the param- 
eters derived at Zpi,oi and those at Zspec- The numbers in the parenthesis 
shows the fraction at which the two parameters agree within Icr. The 
middle panel shows the differences we obtain if we do not use the log- 
ical age constraint. The numbers in the bottom panel mean the median 
differences between parameters from the secondary burst models and 
those from the fiducial t models. 
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Fig. 5. SFR from the SED fits plotted against SFR from the 
MIPS data in GOODS (circles), MIPS in PKSl 138 (squares) for 
galaxies around the radio galaxy redshift (z = 2.15). The double 
circles and squares show the spectroscopically confirmed mem- 
bers in GOODS and PKSl 138, respectively, and the crosses are 
confirmed outliers. Here we plot only galaxies used in the main 
analysis. The dashed line shows SFRsed = SFRmips and the dot- 
ted lines show 1 dex offsets. 



the observed offset is zero is because the model grids are discrete 
and many galaxies remain in the same grid. SFRs do not have 
grids as they are normalized by the observed fluxes of galaxies, 
but most galaxies have the same best-fitting templates and the 
median SFRs remain the same. 
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Fig. 4. Comparisons between the parameters derived at Zsp^c and Zphot- The parameters plotted in the x-axis are derived by fixing 
redshifts at z^pec, and those in the y-axis are derived using redshift as a free parameter The panels show (a) age, (b) star formation 
time scale t, (c) SFR, and (4) dust extinction in Ay. The squares and circles are objects in PKS1138 and GOODS, respectively. The 
thick symbols show the spectroscopic members (i.e., galaxies at z ~ 2.15) and pale symbols show non-members. For clarity, we plot 
the error bars for every 5 objects for non-members. The points are randomly shifted by a small amount to avoid overlapping. 



4.3. SFR from SED Fits 

Next, we check the accuracy of the SFRs from the SED fits. 
There are many ways to derive SFRs of galaxies, but here we 
take those from the warm dust emission probed by the MIPS 
observations. We ta ke the conve r sion fr om the MIPS flux to 
the total IR flux by iReddy et al.1 (l2006h an d derive SFRs us- 
ing the formula given bv iKennicuttI 01998). We have a deep 
MIPS 24yum ima ge in both PK S 1138 (Seymour et al. 2007) and 
GOODS (ISantinietal.ll2009l) . reaching down to ~ 40Mo yr ' 
and ~ 15Mo yr"^ respectively, at z = 2.15. It should be pos- 
sible to measure SFRs from various spectroscopic observations 
carried out in the GOODS field, but the spectra are not always 
public. 

We plot in Fig.|5]SFRs from the SED fits against SFRs from 
MIPS. Note that we plot only galaxies that are detected in MIPS 
and are at Zpimt = 2.15 within 2cr used in the main analyses 
presented in the following sections. For the GOODS sample, 
the agreement between SFRsed and SFRmips is relatively good. 



although there are some outliers. Roughly 70% of the galaxies 
have consistent SFRs from the SED fits and MIPS photometry 
within 1 dex at Icr. Note that SFRmips itself must have a signifi- 
cant error given th e large scatter in mid-far IR SEDs of galaxies 
(lElbaz et al.ll200l . 

The PKS 1 138 sample has a similar rate of catastrophic fail- 
ures to the GOODS sample - 6 out of 24 galaxies (~ 25%) plot- 
ted fall below MIPS SFRs by more than 1 dex, although only 
one of them is inconsistent at Icr. SFRs for the PKS 11 38 galax- 
ies have larger errors compared to the GOODS SFRs. This is 
because some of our optical images in the PKS 1138 are shal- 
low compared to the GOODS images. The optical images probe 
rest-frame UV at z ~ 2 and they are crucial to pin down accurate 
SFRs. 

We have visually inspected the galaxies with SFRsed + 
cr(SFR) < IOMq yr"' (i.e., the upper error is lower than 
IOMq yr"') used in the main analysis on the MIPS images. 75% 
(9 out of 12) of galaxies with low SFRs in GOODS are detected 
in MIPS and have SFRmips ^ 40 Mq yr '. The other three are 
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also detected but have < 20 Mq yr '. This high fraction of the 
MIPS detection illustrates the fundamental difficulty in distin- 
guishing dusty starbursts from passive galaxies. But, this can be 
due to wrong photometric redshifts, AGN contamination, or ob- 
ject blending/confusion in the MIPS image due to the large PSF. 
We ha ve looked at X-ray point sources in the field jLehmer et alj 
l2005h and found that none of the low SFR galaxies are detected. 
We note in passing that 3 out of 50 GOODS galaxies used in the 
main analysis are detected in X-rays, but they do not affect our 
results. 

In contrast to GOODS, the PKS1138 galaxies used in the 
main analysis show a difference between low and high SFR 
galaxies - only ~ 20% of galaxies (4 galaxies out of 21) with 
SFRsED + o-{SFR) < IOM0 are detected in the MIPS image. 
Our MIPS image reaches SFRmips ~ 4OM0 at z = 2.15 and if 
all the low SFR galaxies in PKSl 138 are forming stars at similar 
rates to the GOODS galaxies, they should have been detected. 
We will show that PKS1138 galaxies have lower SFRs than 
GOODS galaxies in the next section and the observed lower fre- 
quency of the MIPS detections in PKS1138 gives a strong sup- 
port to this result. We note that Galametz et al. (2009) suggested 
an increased fraction of AGNs in proto-clusters and some of the 
PKS1138 galaxies with low SFRs may host an AGN, boosting 
the mid-IR flux. In fact, one of the low S FR galaxies that is use d 
in the main analysis is detected in X-rays ( Pentericci et al.l2002h . 
This is the only object detected in X-rays and used in the main 
analysis and this object does not affect our result in any way. 

We admit that the correlation between SFRsed and those 
measured from other methods is not particularly good, and we 
do not attempt to discuss precise SFRs of individual galaxies. 
Instead, we classify galaxies into two coarse classes by their 
SFRs. We define high SFRs as SFRsed > IOMq yr ', and low 
SFRs are lower than that. We do not try to further subdivide the 
SFRs. However, we will later introduce another class - starbursts 
- as S FRsED > IOOM0 yr"' just to illustrate starbursting popu- 
lations at the time of cluster formation at very high redshifts. 

4.4. Secondary burst models 

As we discuss later, cluster galaxies may have experienced many 
mergers in the past and the r models may not reproduce their 
star formation histories very well. We have implemented a sec- 
ondary burst on top of the r models and checked how the results 
change. We assume an instantaneous burst and add a secondary 
burst either at 0.2, 0.5, 1, 2, and 5 Gyr after the onset of star for- 
mation and the burst adds either 10% or 50% of stellar mass of 
the galaxy. 

The bottom line of this exercise is that our conclusions re- 
main unchanged. The bottom panel of Table|2]quantifies the dif- 
ference between the secondary burst models and the t models. 
We observe no strong differences in age and t. The agreement of 
SFRs between the two models is not particularly good, but our 
results remain unchanged because the systematic differences are 
very small. The dust extinction shows a poor agreement and a 
systematic offset in Ay. But, the amount of the offset is similar 
in PKS1138 and GOODS, and our conclusion will therefore re- 
main the same as we discuss only relative differences between 
the two samples. 

The secondary burst models are probably more realistic than 
the simple r models. However, the addition of the secondary 
burst raises the number of free parameters in our SED fitting 
from 5 to 7 (or 8 if we allow the star formation time scale of 
the secondary burst to vary), while we only have 1 1 broad band 
photometry points. As this introduces too many degeneracies be- 



tween the parameters and the secondary burst complicates the 
interpretation of the results, we will consider only the r models 
for our main analysis in the following sections. 



5. Results 

We now move on to present the results from the SED fits. We 
stress again that we use the two samples with almost the same 
wavelength samplings, feed them to the same code with the same 
templates, and discuss the relative differences between the two 
samples. 

In what follows, we study galaxies with K, < 22.5, detected 
in more than 5 bands, and being consistent with Zphoi = 2.15 
within 2cr (note the PKSl 138 radio galaxy is atz = 2.15). We ex- 
clude galaxies with poor fits (xl > 3) as they likely have wrong 
photometric redshifts as discussed in the last section. We take 
the same sample definition in both PKSl 138 and GOODS, so 
that we can make a fair comparison. The physical parameters 
we focus on in this section are age, star formation time scale 
(t), dust extinction, and SFR. We first present evidence that the 
PKSl 138 is indeed an over-density region as suggested by the 
earlier studies. We then move on to discuss physical properties 
of galaxies from the SED fits as a function of environment. 

5.1. Over-density of galaxies in PKS1 138 

Fig. |6] compares the distributions of galaxies around the radio 
galaxy redshift in the PKSl 138 and GOODS fields. We could 
not obtain a good enough photometric redshift for the PKSl 138 
radio galaxy, Zpiwi = '^■^^-o'm (i-^-' i*- ^^ ^-^ '^ away from the 
spectroscopic redshift) and it goes out of our main sample. But, 
the radio galaxy is a powerful AGN and our photometric redshift 
is relatively good for such an extreme object. We show the radio 
galaxy in Fig. |6]just to illustrate the location of the radio galaxy. 

The galaxies around PKSl 138 are arrayed in a filamentary 
structure in the east-west direction, and the radio galaxy is lo- 
cated at the center. There is an over-density of galaxies around 
the central radio galaxy. Interestingly, many of the galaxies 
around the radio galaxy show relatively low SFRs. On the other 
hand, the galaxy distribution in the GOODS field is quite sparse 
and there is no concentration of low SFR galaxies. The surface 
densities of galaxies (Ks < 22.5) are 3.12 + 0.44 arcmin"^ and 
0.64+0.09 arcmin-2 in PKS 1138 and GOODS, respectively, sug- 
gesting an over-density of a factor of ~ 5 in P KSl 138, which 
is in rough agreement with iKurk et all (l2004bi) who found a 
factor of ~ 4 over-density of Ha galaxies. The over-density is 
more pronounced for bright galaxies ; 0.75 + 0.22 arcmin"^ and 
0.08 + 0.03 arcmin-2 for K, < 21.5 (i.e., a factor of 9). AU 
the above numbers for PKSl 138 are averaged over the entire 
probed field and the over-density can be further pronounced if 
we take only the central part of the PKSl 138 field - the den- 
sity of PKSl 138 galaxies within 1' from the radio galaxy is 
10.2 ±1.8 arcmin"^ for K, < 22.5 (i.e., a factor of 16) and 
1.59 + 0.71 arcmin-2 for Ks<2l.5 (i.e., a factor of 20). 

One may wonder that nearby fore-/background structure 
may be contaminating the structure we observe in PKSl 138. 
But, the contamination should be small because the observed 
galaxy distribution is simi lar to the distributi on of Ha emitters at 
the radio galaxy redshift (iKurk et al.ll2004bl) - weak concentra- 
tions of galaxies on the West and South-East of the radio galaxy 
are both seen in the distributions of photo-z selected and Ha se- 
lected galaxies. Also, structures at different redshifts are unlikely 
to make a coherent structure around the radio galaxies with a 
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Fig. 6. Distributions of galaxies with K^ < 22.5 at Zphor = 2.15 within 2cr in the PKS1138 field (left) and GOODS field (right). The 
filled circles show galaxies with low SFRs. To be conservative about red galaxies, we define low SFR galaxies as SFR + cr(SFR) < 
10 Mq yr ' (i.e., the upper error is lower than 10 M© yr '). The open circles are galaxies with high SFR. The stars show galaxies 
fit with AGN templates. The double circles and crosses mean spectroscopically confirmed members and outliers, respectively. The 
sizes of the symbols correlate with the K, band luminosity. The top and right axes show comoving scales at z = 2.15. 



concentration of low SFR galaxies. IZirm et all (l2008l) found a 
concentration of red galaxies selected in J - H colors around 
the very vicinity of the radio galaxy (< 10"), we do not observe 
such a strong concentration of red galaxies, but it might be a con- 
centration of faint galaxies (we use Ks < 22.5, while they used 
H < 24.5, roughly corresponding to Ks < 24). 

Fig. Q plots color magnitude diagrams of galaxies in the 
two fiel ds. We model the location of the red sequence using the 
updated Bruzual & Chariot! (l2003h p opulation synthesis m odels 
following the procedure described in iLidman et al.l (12008 ). The 
most striking trend is that there is a hint of the red sequence, 
which is a ubiquitous feature of galaxy clusters at lower red- 
shifts, at Ks < 21 in PKS1138. The bright red galaxies are 
around the model red sequence formed at Zf ~ 4. This sug- 
gests that these red galaxies formed the bulk of their stars around 
that redshift. The possible red sequence in PKS1138 is popu- 
lated half by galaxies with low SFRs (< 10 Mq yr ') and half 
by galaxies with high SFRs, suggesting that the red sequence is 
being formed. We note that spectroscopically confirmed contam- 
inant galaxies tend to be blue, and thus the possible red sequence 
in PKSl 138 is not due to contamination from fore -/background 
galaxies. In contrast to PKSl 138, there is no clear sign of a red 
sequence in the GOODS field, and there are very few galaxies 
brighter than Ks =21. 

To further illustrate the over-density of PKSl 138, we plot 
stellar mass functions in Fig.|8] We do not derive volume densi- 
ties of galaxies as we are probably looking at a coherent struc- 
ture in redshift in PKSl 138, whose size we do not yet know 
precisely. We simply derive a surface density. The over-density 
of PKSl 138 is evident - PKSl 138 hosts galaxies more densely 
than GOODS. We will focus on galaxies more massive than 



1O"M0 in the next section, for which our samples are nearly 
complete in both PKSl 138 and GOODS. 

To sum up, we observe the significant over-density of galax- 
ies and the forming red sequence in the PKSl 138. This suggests 
that there is a proto-cluster (or perhaps a real cluster) around the 
radio galaxy. But, we have to wait for intensive spectroscopic 
follow-up observations to confirm it. In contrast to PKS 1 1 38, we 
do not observe any possible (proto-)clusters nor red sequence in 
the GOODS field and the galaxy density is lower. Therefore, the 
two samples probe different environments at z ~ 2 and they pro- 
vide us with a unique opportunity to quantify the dependence of 
galaxy properties on density at this high redshift. 

5.2. Physical parameters from the SED fits 

We are now going to look into differences in galaxy properties 
between the two environments. We quantify four fundamental 
properties of galaxies - age (which is time since the onset of star 
formation), star formation time scale, SFR, and dust extinction 
as summarized in Fig. |9] We are aware of the extremely uncer- 
tain nature of various parameters gone into the SED fits. For ex- 
ample, initial mass function (IMF) is one of the most uncertain 
assumptions. A number of IMFs are suggested in the literature, 
each of which produce systematically different results, and to 
make matters worse, IMFs may depend on mode of star forma- 
tion, e.g., IMFs could become top-heavy during starbursts. Also, 
the r models we adopt may not be a good approximation of star 
formation histories of real galaxies, we assume all the stars have 
the same metallicities, and etc. Again, all the uncertainties are 
inherent in the both samples and we discuss only the relative 
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Fig. 7. J - Ks plotted against Kg. The dots are all the detected 
galaxies in each field and large symbols are galaxies around the 
cluster redshift (Zphoi = 2.15 within 2cr). The meanings of the 
symbols are the same as in Fig.|6] Namely, the filled circles and 
open circles are galaxies with low and high SFRs, and the stars 
are AGNs. The red bright galaxy with /Tj = 18 in the left panel 
is the PKS1138 radio galaxy. The vertical dashed lines show 
Ks = 22.5. and the slanted lines show the 3cr limiting colors. 
The dotted lines show the location of the model red sequence 
formed at z = 3, 4 and 5 from bottom to top. 



n PKS1136 
o GOODS 



S?<10 



5x10' 



" 2xl0' 

stellar ["e' 



5x10' 



10' 



"oJ 



Fig. 8. Stellar mass functions of galaxies in PKS1138 (open 
square) and GOODS (open circle). The vertical dashed line 
shows our stellar mass limit of 10"Mq. The error bars show the 
Poisson errors. 



differences between the two samples, which must be the most 
robust results from the SED fits. 

We discuss the four parameters derived from the SED fits. 
Since all the parameters are correlated with stellar masses of 
galaxies and we plot them as a function of stellar mass. Note 
that we miss a fraction of galaxies with masses lower than 
IO^Mq in PKS1138, and we focus on galaxies more massive 
than IO^'Mq. This is a conservative mass cut and our results 
below do not suffer from any strong incompleteness effects. But, 
we keep lower mass galaxies in the plots below to illustrate the 
mass dependency of the four parameters. 

AGE (Fig.lTOh) : The age distributions are not significantly dif- 
ferent. We find that the fractions of galaxies older than 2 Gyr 
are 0.52 + 0.19 and 0.42 + 0.16, respectiveljB Recall that the 



^ As discussed in Section 4.2, we do not use galaxies with poor fits 
(xl > 3). The fractions do not change if we use them; 0.59 ±0.18 
and 0.42 ± 0.15 for PKS 1138 and GOODS, respectively. The fractions 




t [Gyr] 

Fig. 9. The parameters in the models. We obtain star formation 
time scale (r), age, SFR and dust extinction (which is not plot- 
ted) from the fits. We then derive the initial SFR (SFRq) and the 
epoch when galaxies form 80% of their stars (fo g) and discuss 
them in Section 6. 



age is the time since the onset of star formation to the observed 
epoch. These face values should be taken with caution because 
we did not take into account the random errors on the age esti- 
mates. We take a Monte-Carlo approach to quantify the differ- 
ence. We let each data point fluctuate by its error and perform 
the Mann- Whitney test for each realization. We repeat the simu- 
lation for 100,000 times and find that the median probability for 
the two distribution being drawn from the same parent popula- 
tion is 30%. We obtained < 5% probabilities only in 3% of the 
realizations. There is no strong evidence for age differences. 

We note that Steidel et al. (2005) studied a spectroscopically 
confirmed redshift spike at z = 2.3 and found that galaxies in 
the spike are older than those outside the spike. Their result may 
appear inconsistent with ours. But, their sample is a collection 
of star forming galaxies, while we probe the entire galaxy 
population. A fair comparison cannot be made. 

STAR FORMATION TIME SCALE, r (Fig.[lOj}) : PKSl 138 
galaxies tend to have shorter formation time scales. The frac- 
tions of galaxies with t < 0.5 Gyr are 0.70 + 0.23 and 0.43 +0.16 
in PKSl 138 and GOODS, respectively. The Mann- Whitney 
test suggests that the t distributions are likely different (median 
probability of ~ 4%). We obtained < 5% probabilities in 52% 
of the Monte-Carlo realizations. As seen in the GOODS data, 
the fore-/background contaminants do not favor any particular 
T, and the observed difference is not due to the contamination in 
PKSl 138. Combining with the age distribution, most galaxies 
in both fields have age > t, suggesting that they already formed 
the bulk of their stars by the time of the observation. The shorter 
formation time scale of PKSl 138 is one of the most important 
results in this paper, and we will further discuss it in the next 
section. 

SFR (Fig. ITOk) : PKSl 138 galaxies tend to show weaker star 
formation and lower specific SFRs. The fractions of galaxies 
that have SFR < 10 Mq yr"' are 0.74 ± 0.24 and 0.33 ± 0.14 
in PKSl 138 and GOODS, respectively. The Mann- Whitney 
test suggests that they are not from the same parent distribution 
(0.7%). The fore-/background contamination tends to have high 



quoted below for the other parameters also remain essentially the same: 
0.72 + 0.21 and 0.46 ± 0.16 for t <0.5 Gyr, 0.79 ± 0.22 and 0.38 ± 0.14 
for SFR < 10 Mo yr"' , and 0.66 ± 0. 19 and 0.23 ± 0. 10 for Ay < Imag. 
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Fig. 10. Physical parameters from SED fits plotted against stellar mass for PKS1138 and GOODS in the left and right panels in 
each plot, respectively. The vertical dashed lines show the stellar mass limit, below which we are incomplete. The double-circles and 
crosses are spectroscopically confirmed members and fore -/background galaxies, respectively. For objects going outside the plotted 
range, we use the double arrows and open arrows for confirmed members and outliers, (a) top-left: Age plotted against stellar 
mass, (b) top-right: Star formation time scale (t) plotted against stellar mass, (c) bottom-left: SFR plotted against stellar mass. The 
slanted dashed lines show the specific star formation rates defined by SFR / stellar mass, (d) bottom-right: Dust extinction plotted 
against stellar mass. We plot the extinction in units of Ay. 



star formation. They unlikely to contribute to the observed 
difference. We recall that we obtained a lower frequency of 
MIPS detections in PKS1138 than in GOODS, which gives a 
further support to the lower SFRs in PKS1138. 

DUST EXTINCTION (Fig.HSl) : PKS1138 galaxies tend to 
have less dust, as expected from their low SFRs. The fraction 
of PKS1138 galaxies with Ay < 1 is 0.61 + 0.21 (Ay = 1 is a 
typical amount of dust in local star forming galaxies). Most of 
galaxies in GOODS are very dusty, and there see ms a correla- 
tion b etween stellar mass and dust extinction (e.g., Reddv et al. 
l200dl The fraction of Ay < 1 galaxies is 0.17+0.09. The Mann- 
Whitney test shows that the Ay distributions are different (0. 1 %). 

To sum up, we find that galaxies in PKS1138 and GOODS 
have different properties. Interestingly, the trend we observe in 
z < 1 clusters still qualitatively persists even at this high redshift 



— we observe galaxies with lower SFR and less dust in higher 
density regions, where we see clearer red sequence. The envi- 
ronmental dependence of galaxy properties is already in place at 
z = 2, at least partly. We will further extend the discussion and 
quantify the environmental dependence of galaxy formation in 
the next section. 



6. Discussion 

From the extensive SED fits of galaxies in PKS1138 and 
GOODS, we have found that massive galaxies in PKS1138 tend 
to have 

1 . similar ages 

2. shorter star formation time scales 

3. lower star formation rates 

4. lower amounts of dust 
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compared to those in GOODS. The combination of first and sec- 
ond points is interesting. Galaxies in PKS1138 and GOODS 
start forming stars at a similar epoch in a statistical sense, but 
PKS1138 galaxies form more rapidly. We recall that our defi- 
nition of age is the time since the onset of star formation. The 
third and fourth points are basically by-products of the first and 
second points. At z = 2, PKS1138 galaxies have already un- 
dergone intense star formation and their SFRs are rapidly de- 
clining, while GOODS galaxies are still actively forming stars 
due to their longer star formation time scales. This results in 
the lower SFRs in PKS1138 (third point). The fourth point can 
then be easily understood given the correlation between between 
SFR and dust amou nt (lower SFR galaxies have less dust, e.g., 
iHopkins et alll2003h . 

Cluster galaxies have a shorter star formation time scale - 
this is the same tren d as observed in z < 1 clusters. For example, 
iGobat et al] (l2008l) showed the same trend in a z - 1.2 clus- 
ter based on a photo-spectroscopic analysis. It is striking that 
the trend holds even at z = 2.15. The difference in the star for- 
mation time scale suggests that cluster and field galaxies may 
form in different ways. Let us discuss the formation and evolu- 
tion of cluster and field galaxies in detail. We first introduce two 
parameters to further quantify the galaxy formation. We then ex- 
tend the discussion and address the origin of the environmen- 
tal dependence of galaxy properties observed at lower redshifts. 
We note a caveat here that our results are based on one (proto- 
)cluster only and the trends we observe may not represent global 
trends at z = 2. A larger sample of z = 2 clusters should be 
investigated to draw a global picture. 



6.1. The environmental dependence of galaxy formation 

We have age, star formation time scale, and SFR for each galaxy. 
Assuming the exponentially decaying star formation histories, 
we can estimate the SFR at the onset of star formation, which 
we call initial SFR (SFRq; see Fig. [Hi. Fig.Oplots SFRq as a 
function of stellar mass. The cluster galaxies tend to have higher 
SFRo- the fractions of SFRo > lOOOMoyr' are 0.61 ±0.21, and 
0.39 + 0.15 in PKS1138 and GOODS, respectively. The Mann- 
Whitney test supports this difference (the median probability is 
2% and we obtain < 5% probability in 74% of the realizations). 
PKS 1138 galaxies have experienced more intense formation his- 
tories. Of course, all this is based on the assumption of the expo- 
nentially declining SFRs, and we have ignored effects of galaxy- 
galaxy mergers, which we will discuss later. 

If galaxies continue to form stars following the exponential 
decay, we can derive the time when the galaxies form bulk of 
their stars. We plot in Fig. [12] the epoch when galaxies form 
80% of the stars that they would have at z = (fo.s; see Fig. 
|9]l. PKS 1138 galaxies typically formed around z ~ 3 or higher, 
while GOODS galaxies typically formed below z ~ 3. We recall 
that we obtained z/ ~ 4 from the location of the red sequence 
in PKSl 138 in Section 5.1, which is in agreement with what we 
find here. We find that 0.61 + 0.21 of galaxies in PKSl 138 form 
the bulk of their stars at z > 3, while the fraction is 0.17 ± 0.09 
in GOODS. The Mann- Whitney test suggests that they are likely 
different - the median probability is ~ 8 % and we obtained < 5 % 
probabilities in 38% of the realization. Interestingly, the ages of 
the cluster and field galaxies are not very different (Fig. [TOk). 
The difference in fo.8 is therefore due to their shorter formation 
time scales. The formation of cluster galaxies is a more intense 
event, they form in a shorter time scale, and they assemble the 
bulk of their stars earlier by ~ 1 - 2 Gyr than the field galaxies. 
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Fig. 13. The median star formation histories of the PKS 1138 and 
GOODS galaxies. The shaded area shows 68 percentile intervals. 



To summarize, we measure the median r and SFR for > 
1O"M0 galaxies and illustrate the differences in star forma- 
tion histories of galaxies in PKSl 138 and GOODS in Fig. [13] 
This will the summary plot of the paper As discussed above, 
PKSl 138 galaxies have higher SFRo and shorter t. PKSl 138 
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experience much more intense galaxy formation at early times 
than GOODS galaxies. The galaxy assemblies are completed in 
a short time, and by the time we observe them (~ 2 Gyr after the 
onset of star formation as shown in Fig.fTOh) their SFRs become 
low. The plot summarizes the points listed at the beginning of 
this section (except for the fourth one, which is a result of SFR- 
dust correlation). 

The plot indicates that PKS1138 galaxies and GOODS 
galaxies did not form in the same way. The formation of cluster 
galaxies is a more intense event. But, how can massive galaxies 
form on such a short time scale? In general, star formation ac- 
tivities have negative feedback to themselves. Hot, young stars 
ionize surrounding gas, which prevents further star formation. 
Also, supernova explosions of massive stars give negative feed- 
back. If galaxies undergo an AGN phase, it might also give nega- 
tive feedback. The formation of massive galaxies on a very short 
time scale may be difficult to understand in this respect, and that 
makes us speculate that these massive galaxies might have taken 
another route - mergers. 

Galaxy formation takes place in high peaks of density fluc- 
tuations in the Universe. Clusters that we observe today were 
regions of large-scale over-densities, where many density peaks 
were embedded. On the other hand, field galaxies form in more 
isolated peaks. The current galaxy formation theories predict 
that small galaxies form first and they form progressively mas- 
sive galaxies via mergers. In the early Universe, cluster galaxies 
must have had a lot better chances to merge with other galax- 
ies than field galaxies simply because there are more galaxies 
around them. This accelerated frequency of mergers might be 
the cause of the rapid formation of cluster galaxies. 

A simple scenario of an equal mass merger can be helpful 
here to show that mergers tend to make SFRo higher and push 
?o.8 to higher redshifts. Assuming that two equal mass galaxies 
merge into a single galaxy, the merger event will double the stel- 
lar mass of the galaxy. If the galaxy experiences starburts trig- 
gered by the merger, its stellar mass will be more than double 
and its SFR gets lower after the burst. We then observe this mas- 
sive galaxy with low SFR, and depending on the time elapsed 
after the burst, we would fit a short t to this galaxy in order to 
reproduce its low SFR and large stellar mass. In fact, we find that 
the secondary burst discussed in Section 4.4 reduces the median 
SFRo of the PKSl 138 galaxis by ~ 30%, although it is stifl high, 
SFRq = I6OOM0 yr"', which might indicate that galaxies expe- 
rienced more than one merger. On the other hand, the S FRq of 
the GOODS galaxies remain the same within 10%. 

We are not yet sure if PKS 1 1 38 is a collapsed cluster or cur- 
rently collapsing proto-cluster. But, even if it is a cluster, it must 
be a young system, and it is unlikely that galaxies had enough 
time to be affected by 'nurture' effects such as ram-pressure 
stripping. We may suggest that mergers in the early times of the 
cluster formation may be one of the key processes to establish 
the environmental dependence of galaxy properties observed at 
lower redshifts. We will further pursue this point in the following 
section. 



6.2. The Nature effects in the environmental dependence of 
gaiaxy properties 

Galaxy properties such as colors and morphology are known to 
depend on environment in which galaxies reside. This environ- 
mental dependence is shaped by two effects - nature and nurture 
effects. That is, both how galaxies form and how they evolve are 
important. 



Fig. \T2\ shows that only ~ 1 Gyr has passed since the 
PKSl 138 galaxies formed bulk of their stars to the observed 
epoch. Assuming that PKS1138 is a virialized sy stem, we ap- 
ply a 2cr-clipped gapper method ( Beers etal.ll 1 990h and obtain a 
velocity dispersion of 400 km s ' using spectroscopic redshifts 
within 2 arcmin, co rresponding to a physi c al scale of 1 Mpc, 
from the literature ('Penter icci eTaP 120001: iKurk et al.1 l2004bt 
Croft et al. 2005). The virial radius (^200) of the cluster is 0.32 
Mpc, giving a crossing time scale of ~ 1 Gyr. Therefore, the 
time elapsed since the bulk formation of the PKS 11 38 galaxies 
is comparable to the crossing time scale of this cluster. This sug- 
gests that nurture effects may not have had enough time to fully 
work. They may have affected a fraction of galaxies, but they 
probably could not change the average properties of galaxies. 
We may be witnessing the nature effects in shaping the environ- 
mental dependence in PKS 11 38. 

We admit that it is not very straightforward to classify nature 
and nurture effects at this high redshift. Now, let us introduce two 
kinds of mergers; early-epoch mergers and late-epoch mergers. 
We refer to mergers occur during the first collapse of clusters as 
early-epoch mergers, and those occur afterwards as late-epoch 
mergers. This is to sort out the two effects, early-epoch merg- 
ers being nature effects (initial conditions) and late-epoch merg- 
ers being nurture effects (environmental effects). Of course, nur- 
ture effects also include ram-pressure stripping, harassment, etc. 
On the other hand, nature effects do not include these processes 
driven by the deep potential well or intracluster medium. 

In this classification, the environmental dependence of 
PKSl 138 is likely due to early-epoch mergers because ~ 1 Gyr 
time is probably not enough for nurture effects to fully work. 
Galaxy clusters at z ~ 1 do not give useful information about 
nature effects because nurture effects have had enough time to 
fully operate (~ 4 Gyr since the formation epoch) and it is not 
straightforward to disentangle the two effects. Our results sug- 
gest that nature effects are strong and they form the basis of the 
environmental dependence of galaxy properties. A way to probe 
the significance of the nature effects is to quantify morphology 
of the PKSl 138 galaxies. If early-epoch mergers are an impor- 
tant effect, then we expect to observe more signs of recent in- 
teractions in PKSl 138 than in GOODS (see below). Also, we 
expect to observe post-starburst galaxies if galaxies undergone 
interaction triggered starbursts. A near-IR spectroscopic follow- 
up campaign of PKSl 138 is currently underway. We will be able 
to study spectral properties of the PKSl 138 galaxies, which al- 
low us to look deeper into star formation histories. 

To sum up, the strong nature effects may shape the environ- 
mental dependence of galaxy properties. An ultimate goal of en- 
vironment studies will be to quantify the relative contribution of 
nature and nurture effects. But, that will require statistical work 
on z ~ 2 (proto-)clusters. Although the possible significance 
of nature effects we suggest here is based only on one (proto- 
)cluster PKSl 138, let us further discuss it. It has an interesting 
implication for the build-up of the cluster red sequence. 



6.3. The massive end of the duster red sequence 

The sequence of red early-type galaxies is a ubiquitous fea- 
ture of galaxy clusters. Over the last few years, there is an ac- 
cumulating amount of evidence in the literature that the clus- 
ter red sequen ce grows from the mass i ve end to the low- 
mass end (e.g., Tanaka et al. 2005, '2007"; 'Kov ama et al.l 120071 : 
Tanaka et al. 2008; Gilbank et al. 2008). The massive end of the 
red sequence cannot be formed via a simple fading of blue, star 
forming galaxies because such massive blue galaxies do not ex- 
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ist even at z ~ 2 (see Fig.|2]and discussions in lFaber et alfcOOTh . 
We need mergers to form it. 

An interesting point here is that we observe the brightest 
tip of the red seq uence in PKS1138, which was also noted by 
IZirm et aTl (|2008|) . The massive end of the red sequence in a 
young system - this might be due to early-epoch mergers. Early- 
epoch mergers might have formed very massive galaxies during 
the first gravitational collapse of clusters. Early-epoch mergers 
should occur more frequently in cluster environments than in the 
field, and that helps explain why we do not observe red sequence 
in GOODS . As suggested by Zirm et al. ( 2008 ), the red sequence 
in PKS1138 may be being formed or just formed at the time of 
observation. This view is further supported by our observation 
that roughly half of the red sequence galaxies have high SFRs 
(Fig. [Tot . This formation redshift of z ~ 2 is in line with pre- 
dictions from the build-up of the red sequence observed in lower 
redshift clusters (iTanaka et al.ll200'7l) . 

We present in Fig. [14] ACS /-band images of the bright red 
galaxies in PKS1138 {K, < 21 and J - Ks > 1). Here we only 
briefly discuss their morphologies and a detailed study will be 
presented elsewhere (Zirm et al. in prep). Interestingly, a half 
of the galaxies show disturbed morphology and/or have nearby 
companions, lending a support to the picture of the accelerated 
mergers in clusters. Two out of three apparently disturbed galax- 
ies are detected in MIPS. Early-epoch mergers during the first 
collapse of clusters may form the brightest end of the red se- 
quence, and at the same time, form the basis of the environmen- 
tal dependence of galaxy properties. 

The low-mass end of the red sequence is built up at later 
times. Nurture eff'ects may come in there. Once clusters form, in- 
tracluster gas and deep potential field, in addition to late-epoch 
mergers, can affect galaxies and terminate their star formation 
activities. In this way, moderate-low mass galaxies could be- 
come red and form the low-mass end of the red sequence. The 
combination of nature and nurture eff'ects may have conspired 
to produce the down-sizing behavior of the observed build-up of 
the red sequence (e.g.,iTanaka et al.i ii2005i) . 

6.4. Proto-clusters at higher redshifts 

Finally, we finish the discussion with a forecast for future (proto- 
)cluster studies. We can go back in time and see how the 
PKS 1138 field looked like at higher redshifts as we have star for- 
mation histories of individual galaxies. A caveat of course is that 
we ignore all the early-epoch mergers occurred before the time 
of the observation (the only way to recover that information is to 
resolve the galaxies into individual stars). Another caveat is that 
we cannot track spatial positions of galaxies back in time. But, 
PKS 1 138 is likely a collapsing/collapsed system and it may well 
have been an over-density region already at z = 4. 

Fig. [15] shows distribution of galaxies and their star forma- 
tion rates evolved back to z = 4. Compared to Fig. [6] galaxies are 
more actively forming stars on average. Interestingly, galaxies 
in PKS 1 138 are more actively forming stars and there are more 
starbursting galaxies with SFR > 100 Mq yr ' than in GOODS. 
This is in stark contrast to Fig. [6] where we saw that PKS 1 138 
galaxies have lower SFRs than those in GOODS. We still ob- 
serve a hint of a galaxy over density in PKS 1 138 at z = 4. The 
plots suggest that, as we approach the formation epoch of clus- 
ters (i.e., early phase of the gravitational collapse to a massive 
cluster halo), we expect to observe an over density of starburst- 
ing galaxies. Along with starbursting galaxies, low SFR galaxies 
already appear in PKS 1138, while such galaxies are extremely 



rare in GOODS. At this redshift, low SFR galaxies and starburst- 
ing galaxies may co-exist in a forming cluster 

It is not easy to show how early-epoch mergers change 
the picture we see here. But, we still expect to observe star- 
bursting galaxies triggered by early-epoch mergers in collapsing 
clusters with a higher over-densities of lower-mass, pre-merger 
galaxies. We deem that observations with existing/future sub- 
millimeter arrays would be able to discover many forming clus- 
ters at high redshifts. In fact, some submm observations of dis- 
tant radio galaxies have found possible over-densities of dusty 
starburst p opulations around them (e.g.. iDe Breuck et al.ll2004l: 
iGreve et al.ll20 07). Full wavelength observations will be essen- 
tial for future proto-cluster studies at very high redshifts. 



7. Summary 

We have studied the environmental dependence of galaxy prop- 
erties at z ~ 2 based on the multi-band data available in the 
(proto-)cluster field PKS1138 and in GOODS. We have per- 
formed the extensive SED fits with a special care for systematic 
biases between the two samples. The results from the careful 
SED fits suggest that the environmental dependence is at least 
partly in place at this high redshift. 

We have first shown that PKS 1 138 is indeed an over-density 
region with an excess of red galaxies, forming the brightest tip of 
the red sequence, compared to GOODS. The red galaxies tend to 
cluster around the radio galaxy. These results support the claim 
that PKS1138 is a (proto-)cluster at z = 2.15. Interestingly, the 
red sequence is populated both by low SFR galaxies and high 
SFR galaxies, suggesting that the red sequence is being formed. 

We then have looked into detailed properties of galaxies de- 
rived from the SED fits. PKS 1 138 galaxies have similar age (we 
define age as time since the onset of star formation), shorter 
star formation time scale, lower SFR, and less dust compared 
to those in GOODS at similar redshifts. The averaged star for- 
mation history of the PKS 1138 suggests that the cluster galax- 
ies form on a shorter time scale and they form the bulk of their 
stars ~ 1 Gyr earlier than the field galaxi es, which is consis- 
tent with lower re dshift observations (e.g., [Thomas et al.ll2005l; 
lGobatetal.ll2008l) . 

The environmental dependence of galaxy properties should 
be shaped both by nature and nurture effects. The observed en- 
vironmental dependence at z = 2.15 suggests that nature effects 
may be a strong effect as the PKS 1138 (proto-)cluster is likely 
a young system. Possibly an accelerated rate of mergers in col- 
lapsing clusters may have played a role there and they may also 
be the primary cause of the short formation time scale of cluster 
galaxies. However, further studies on z > 2 (proto-)clusters will 
be needed to confirm this picture. We expect that star bursting 
galaxies populate in very high redshift proto-clusters and sub- 
millimeter observations may be useful to confirm very high-z 
systems. 
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Fig. 14. ACS /-band images of the bright red galaxies. 

A R.A. [Mpc] 
10 -10 



10 



A R.A. [Mpc] 




10 



LO 



CJ 



Q 
<] 



o - 



m 









_PKS1138 (evolved back to z=4.0) 


GOODS (evolved back to 

o 


z = 4.0) 


" * <lo8Ma o lO^-^X o 101°-"Mq O>10"Mg 




- 


□ n n 


, 


. ,o 

■ ■ D ■■- . ■■■■■■.■.- 
. .Bo 




':;■',." □. - 


■ fioEffiO 
o 0-- /,■■-;. . 




■o 


. - ^^ 






. SFR<10Mq yr"^ 




Q 


o lOMg yr"^<SFR<100MQ yr"^ 


■■□... ■. . '■ .. "'.■■- 


.■"■■;'/^' -. : , .- 


D SFR>1OOM0 yr"^ 







o 



u 



- o 



o 

Q 



o 







A R.A. [arcmin 



-5 



A R.A. [arcmin] 



Fig. 15. Same as Fig. |6] but galaxies are evolved back in time to z = 4 assuming the star formation histories obtained from the SED 
fits. The sizes of the symbols correlate with stellar mass. The top and right axes show comoving scales at z = 4. 

obtained from the Hubble Legacy Archive, which is a collaboration between Berlin, E. & Amouts, S. 1996, A&AS, 117, 393 

the Space Telescope Science Institute (STScI/NASA), the Space Telescope Blakeslee, J. P., Franx, M., Postman, M., et al. 2003, ApJ, 596, L143 

European Coordinating Facility (ST-ECF/ESA) and the Canadian Astronomy Bruzual, G. & Chariot, S. 2003, MNRAS, 344, 1000 

Data Centre (CADC/NRC/CSA). Chabrier, G. 2003, PASP, 115, 763 

Chariot, S. & Fall, S. M. 2000, ApJ, 539, 718 

Cooper, M. C, Newman, J. A., Weiner, B. J., et al. 2008, MNRAS, 383, 1058 

Croft, S., Kurk, J., van Breugel, W., et al. 2005, AJ, 130, 867 

References Cucclati, O., lovlno, a., Marinonl, C, et al. 2006, A&A, 458, 39 

n T r- ci T^ ;, n uv, ^, T^ , oon at u,n qt De Breuck, C, Bertoldi, R, CariUi, C, et al. 2004, A&A, 424, 1 

Beers, T C, Flynn, K., & Gebhardt, K. 1990, AJ, 100, 32 



14 



Tanaka et al.: The Environmental Dependence of Galaxy Properties at z = 2 



Doherty et al. 2009, A&A, submitted 

Elbaz, D., Cesarsky, C. J., Chanial, P., et al. 2002, A&A, 384, 848 

Elbaz, D., Daddi, E., Le Borgne, D., et al. 2007, A&A, 468, 33 

Faber, S. M., Willmer, C. N. A., Wolf, C, et al. 2007, ApJ, 665, 265 

Furusawa, H., Shimasaku, K., Doi, M., & Okamura, S. 2000, ApJ, 534, 624 

Galametz, A., Stem, D., Eisenhardt, P. R. M., et al. 2009, ApJ, 694, 1309 

Gebhardt, K., Faber, S. M., Koo, D. C, et al. 2003, ApJ, 597, 239 

Gilbank, D. G., Yee, H. K. C, Ellingson, E., et al. 2008, ApJ, 673, 742 

Gobat, R., Rosati, R, Strazzullo, V., et al. 2008, A&A, 488, 853 

Grazian, A., Fontana, A., de Santis, C, et al. 2006, A&A, 449, 951 

Greve, T. R., Stern, D., Ivison, R. J., et al. 2007, MNRAS, 382, 48 

Gunn, J. E. & Stryker, L. L. 1983, ApJS, 52, 121 

Hatch, N. A., Overzier, R. A., Kurk, J. D., et al. 2009, MNRAS, 395, 1 14 

Hatch, N. A., Overzier, R. A., Rottgering, H. J. A., Kurk, J. D., & Miley, G. K. 

2008, MNRAS, 383, 931 
Hopkins, A. M., Miller, C. J., Nichol, R. C, et al. 2003, ApJ, 599, 971 
Kennicutt, Jr, R. C. 1998, ARA&A, 36, 189 

Kodama, T., Tanaka, 1., Kajisawa, M., et al. 2007, MNRAS, 377, 1717 
Koyama, Y., Kodama, T, Tanaka, M., Shimasaku, K., & Okamura, S. 2007, 

MNRAS, 382, 1719 
Kriek, M., van Dokkum, P G., Franx, M., et al. 2008, ApJ, 677, 219 
Kuntschner, H., Smith, R. J., CoUess, M., et al. 2002, MNRAS, 337, 172 
Kurk, J. D., Pentericci, L., Overzier, R. A., Rottgering, H. J. A., & Miley, G. K. 

2004a, A&A, 428, 817 
Kurk, J. D., Pentericci, L., Rottgering, H. J. A., & Miley, G. K. 2004b, A&A, 

428, 793 
Kurk, J. D., Rottgering, H. J. A., Pentericci, L., et al. 2000, A&A, 358, LI 
Lehmer, B. D., Brandt, W. N., Alexander, D. M., et al. 2005, ApJS, 161, 21 
Lidman, C, Rosati, R, Tanaka, M., et al. 2008, A&A, 489, 981 
Madau, P 1995, ApJ, 441, 18 

Mei, S., Holden, B. R, Blakeslee, J. P, et al. 2009, ApJ, 690, 42 
Miley, G. & De Breuck, C. 2008, A&A Rev., 15, 67 
Miley, G. K., Overzier, R. A., Zirm, A. W., et al. 2006, ApJ, 650, L29 
Nakata, F., Kodama, T, Shimasaku, K., et al. 2005, MNRAS, 357, 1357 
Pentericci, L., Fan, X., Rix, H., et al. 2002, AJ, 123, 2151 
Pentericci, L., Kurk, J. D., Rottgering, H. J. A., et al. 2000, A&A, 361, L25 
Pentericci, L., Roettgering, H. J. A., Miley, G. K., Carilli, C. L., & McCarthy, P. 

1997, A&A, 326, 580 
Pentericci, L., Roettgering, H. J. A., Miley, G. K., et al. 1998, ApJ, 504, 139 
Polletta, M., Tajer, M., Maraschi, L., et al. 2007, ApJ, 663, 81 
Postman, M., Franx, M., Cross, N. J. G., et al. 2005, ApJ, 623, 721 
Press, W. H. & Schechter, P 1974, ApJ, 187, 425 
Reddy, N. A., Steidel, C. C, Fadda, D., et al. 2006, ApJ, 644, 792 
Salpeter, E. E. 1955, ApJ, 121, 161 

Santini, P., Fontana, A., Grazian, A., et al. 2009, ArXiv e-prints 
Schlegel, D. J., Finkbeiner, D. P, & Davis, M. 1998, ApJ, 500, 525 
Seymour, N., Stern, D., De Breuck, C, et al. 2007, ApJS, 171, 353 
Sirianni, M., Jee, M. J., Benitez, N., et al. 2005, PASR 117, 1049 
Springel, V., White, S. D. M., Jenkins, A., et al. 2005, Nature, 435, 629 
Stanford, S. A., Romer, A. K., Sabirli, K., et al. 2006, ApJ, 646, L13 
Steidel, C. C, Adelberger, K. L., Shapley, A. E., et al. 2005, ApJ, 626, 44 
Tanaka, M., Finoguenov, A., Kodama, T, et al. 2008, A&A, 489, 571 
Tanaka, M., Kodama, T, Arimoto, N., et al. 2005, MNRAS, 362, 268 
Tanaka, M., Kodama, T, Kajisawa, M., et al. 2007, MNRAS, 377, 1206 
Thomas, D., Maraston, C, Bender, R., & Mendes de Oliveira, C. 2005, ApJ, 

621, 673 
Venemans, B. P, Rottgering, H. J. A., Miley, G. K., et al. 2007, A&A, 461, 823 
Zirm, A. W., Stanford, S. A., Postman, M., et al. 2008, ApJ, 680, 224 



15 



